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Abstract 

We report an experimental investigation on the influence of an external magnetic field on forced 
3D turbulence of liquid gallium in a closed vessel. We observe an exponential damping of the 
turbulent velocity fluctuations as a function of the interaction parameter N (ratio of Lorentz force 
over inertial terms of the Navier-Stokes equation). The flow structures develop some anisotropy 
but do not become bidimensional. From a dynamical viewpoint, the damping first occurs homo- 
geneously over the whole spectrum of frequencies. For larger values of N, a very strong additional 
damping occurs at the highest frequencies. However, the injected mechanical power remains inde- 
pendent of the applied magnetic field. The simultaneous measurement of induced magnetic field 
and electrical potential differences shows a very weak correlation between magnetic field and ve- 
locity fluctuations. The observed reduction of the fluctuations is in agreement with a previously 
proposed mechanism for the saturation of turbulent dynamos and with the order of magnitude of 
the Von Karman Sodium dynamo magnetic field. 
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Situations were a magnetic field interacts with a turbulent flow are found in various 
domains of physics, including molten metals processing, laboratory flows, and astrophysics. 
The motion of electrically conducting fluid in a magnetic field induces electrical currents, 
which in turn react on the flow through the Lorentz force. The power injected in the flow 
is thus shared between two dissipative mechanisms: viscous friction and ohmic dissipation 
of the induced currents. On the one hand, the situation where the flow is laminar is very 
well understood, and the geometry of the velocity field and of the induced currents can 
be computed analytically. On the other hand, several questions remain open in the fully 
turbulent situation: when a statistically steady state is reached, is the mean injected power 
higher or lower than in the nonmagnetic case? What controls the ratio of ohmic to viscous 
dissipation? How is the turbulent cascade affected by the magnetic field? To adress some 
of these questions, we have designed an experimental device which allows to apply a strong 
magnetic field on a fully turbulent flow. 



When an electrically conducting fluid is set into motion, a magnetic Reynolds number 
Rm can be defined as the ratio of the ohmic diffusive time to the eddy-turnover time. This 
number reaches huge values in galactic flows, but can hardly exceed one in a laboratory 
experiment. The common liquids of high electrical conductivity (gallium, mercury, sodium) 
have a very low kinematic viscosity: their magnetic Prandtl number Pm (ratio of the kine- 
matic viscosity over the magnetic diffusivity) is less than 10 -5 . A flow with Rm of order 
one is turbulent and thus requires a high power input to be driven. For this reason most 
experimental studies have been restricted to low Rm. They were conducted mostly in chan- 
nel flows and grid generated turbulence 

a a. 

The most general observation is that the 
application of a strong magnetic field leads to a steeper decay of the power spectra of the 
turbulent velocity fluctuations at large wavenumbers: the decay goes from a classical k~ 5 ^ 3 
scaling without magnetic field to k~ 3 or k~ 4 for the highest applied fields. In the meantime, 
some anisotropy is developed leading to larger characteristic scales along the applied mag- 
netic field. The phenomenology of these transformations is quite well understood in terms 
of the anisotropy of the ohmic dissipation [3]. One difficulty arises from the fact that the 
boundary conditions can have a strong effect on the turbulence level. In channel flows the 
choice between conducting and insulating walls strongly impacts the flow. A strong external 
magnetic field perpendicular to the boundaries leads to an increase of the turbulence level 
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due to modifications of the boundary layers pj. On the contrary, Alemany et al. observed in 
grid generated decaying turbulence an enhanced decay of the turbulent fluctuations. As far 
as forced turbulence is concerned, Sisan et al. studied the influence of a magnetic field on a 
flow of liquid sodium inside a sphere 4|. However, the flow again goes through a variety of 
instabilities which prevents a study solely focussed on the impact of the magnetic field on 
the turbulent fluctuations: the geometry of the mean flow and of the induced magnetic field 
keeps changing as the magnetic field is increased. 

Several numerical simulations of this issue have been conducted Despite 
the rather low spatial resolution available for this problem, these works show the same phe- 
nomenology: development of anisotropy, steepening of the spectra and trend to bidimension- 
alization of the flow as the magnetic field increases. The direct numerical simulations (DNS) 
also allow to compute the angular flux of energy from the energy-containing Fourier modes 
(more or less orthogonal to the applied field) to the modes that are preferentially damped 
by ohmic dissipation. Once again, the magnetic field can affect the large scale structure of 
the velocity field: the interaction between forcing and magnetic field in the DNS of Zikanov 
& Thess 5, (| leads to an intermittent behavior between phases of roughly isotropic flow 
and phases of bidimensional columnar vortices that eventually get unstable. Vorobev et 
al.Q and Burratini et al.Q] performed DNS in a forced regime which are closely related to 
our experiment. However, the maximum kinetic Reynolds number reached in their studies 
remains orders of magnitude below what can be achieved in the laboratory. An experimental 
investigation at very high kinetic Reynolds number thus remains necessary to characterize 
the effect of a strong magnetic field on a fully developped turbulent cascade. 

The flow under study in our experiment resembles the Von Karman geometry: counter- 
rotating propellers force a flow inside a cylindrical tank. The large scale flow is known to 
have a strong shear layer in the equatorial plane where one observes the maximum level of 
turbulence [9(] . We do not observe any bifurcation of the flow as the magnetic field increases 
nor any drastic change of the large scale recirculation imposed by the constant forcing. 
These ingredients strongly differ from the previously reported experiments and allow to 
study the influence of an applied magnetic field on the turbulent cascade in a given flow 
geometry. This issue arises in the framework of turbulent laboratory dynamo studies such 
as the Von Karman Sodium (VKS) experiment n]] which flow's geometry is similar to ours. 
In the turbulent dynamo problem, one issue is to understand the precise mechanism for the 
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FIG. 1: Sketch of the experimental setup. The curved arrows represent the average large-scale 
motion of the fluid. 

saturation of the magnetic field. Once the dynamo magnetic field gets strong enough, it 
reacts on the flow through the Lorentz force. This usually reduces the ability of the flow to 
sustain dynamo action. An equilibrium can be reached, so that the magnetic field saturates. 
This backreaction changes the properties of the bulk turbulence. In a situation where an a 
effect takes part in the generation of the dynamo, i.e. if the turbulent fluctuations have a 
mean field effect, then the changes in the statistics of turbulence should be involved in the 
saturation mechanism of the dynamo. 

I. DESCRIPTION OF THE EXPERIMENT 
A. The flow 



Our turbulent flow resembles the Von Karman geometry, which is widely used in exper- 
iments on turbulence and magnetohydrodynamics (0, [lo| for example). 8 liters of liquid 
gallium are contained in a closed vertical cylinder of diameter 20 cm and height 24 cm. The 
thickness of the cylindrical wall is 7.5 mm, and that of the top and bottom walls is 12 mm 
(Fig. [T]). This tank is made of stainless steel. Gallium melts at about 30°C. Its density is 
p = 6090 kg/m 3 and its kinematic viscosity is v = 3.11 10~ 7 m 2 /s. Its electrical conductivity 
at our operating temperature (45°C) is a = 3.9 10 6 fi~ 1 m~ 1 . The main difference with the 
usual Von Karman setup is the use of two propellers to drive the flow, rather than impellers 
or disks. The two propellers are coaxial with the cylinder. They are made of 4 blades 
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inclined 45 degrees to the axis. The two propellers are counter-rotating and the blades are 
such that both propellers are pumping the fluid from the center of the tank towards its end 
faces. The large scale flow is similar to the traditional counter-rotating geometry of the Von 
Karman setup. The propellers are 7 cm in diameter and their rotation rate f rot ranges from 
3 Hz to 30 Hz. They are entrained at a constant frequency by DC motors, which drives a 
large scale flow consisting in two parts : first the fluid is pumped from the center of the 
vessel on the axis of the cylinder by the propellers. It loops back on the periphery of the 
vessel and comes back to the center in the vicinity of the equatorial plane. In addition to 
this poloidal recirculation, the fluid is entrained in rotation by the propellers. Differential 
rotation of the fluid is generated by the counter rotation of the propellers, which induces a 
strong shear layer in the equatorial plane where the various probes are positioned. The shear 
layer generates a very high level of turbulence [3|: the velocity fluctuations are typically of 
the same order of magnitude as the large scale circulation. 

A stationary magnetic field Bo is imposed by a solenoid. The solenoid is coaxial with 
the propellers and the cylinder. Bo is mainly along the axis of the cylinder with small 
perpendicular components due to the finite size of the solenoid (its length is 32 cm and its 
inner diameter is 27 cm). The maximum magnetic field imposed at the center of the tank 



B. Governing equations and dimensionless parameters 

The fluid is incompressible and its dynamics is governed by the Navier-Stokes equations: 



where v is the flow velocity, p is the pressure, B is the magnetic field and j is the electrical 
current density. The last term in the r.h.s. of (CD) is the Lorentz force. 

A velocity scale can be defined using the velocity of the tip of the blades. The radius of 
the propeller being R = 3.5 cm, the velocity scale is 2nRf rot . A typical length scale is the 
radius of the cylinder L = 10 cm. The kinetic Reynolds number is then 
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At the maximum speed reported here it reaches 2 10 6 . The flow is then highly turbulent and 
remains so even for a rotation rate ten times smaller. 

In the approximations of magnetohydrodynamics, the temporal evolution of the magnetic 
field follows the induction equation: 

<9B 

_ = V x (v x B) + i]AB . (4) 

C L 

Here rj = (/i o") _1 is the magnetic diffusivity, /Xq being the magnetic permeability of vacuum 



(77 = 0.20 m 2 s _1 for gallium so that Pm = u/rj = 1.6 10~ 6 ) The first term in the 

r.h.s describes the advection and the induction processes. The second term is diffusive and 
accounts for ohmic dissipation. A magnetic Reynolds number can be defined as the ratio of 
the former over the latter: 

Rm = n a27rRLf rot , (5) 
t is a measure of the strength of the induction processes compared to the ohmic dissipation 



11]. For Rm — > the magnetic field obeys a diffusion equation. For Rm ^ 1, it is 
transported and stretched by the flow. This can lead to dynamo action, i.e. spontaneous 
generation of a magnetic field sustained by a transfer of kinetic energy from the flow to 
magnetic energy. In our study, Rm is about 3 at the maximum speed reported here (f rot = 
30 Hz): Induction processes are present but do not dominate over diffusion. As f rot goes 
from 3 to 30 Hz, the magnetic Reynolds number is of order 1 in all cases. 

In our range of magnetic Reynolds number, the induced magnetic field b = B — B is 
weak relative to B (of the order of 1%). The Lorentz force is j x B where j is the current 
induced by the motion of the liquid metal (no exterior current is applied). From Ohm's 
law, this induced current is j ~ o~v x Bo and is thus of order 2iiaRf ro tBQ. One can define 
an interaction parameter N that estimates the strength of the Lorentz force relative to the 
advection term in the Navier-Stokes equation: 

N = " LB o (6) 

Because of the l/f ro t factor, this quantity could reach very high values for low speeds, 
but the flow would not be turbulent anymore. In order to remain in a turbulent regime, 
the smallest rotation rate reported here is 3 Hz. The maximum interaction parameter 
(built with the smallest rotation frequency and the strongest applied magnetic field) is then 
approximately 2.5. Many studies on magnetohydrodynamics use the Hartmann number to 



quantify the amplitude and influence of the magnetic field. This number measures the ratio 
of the Lorentz force over the viscous one. It is linked to the interaction parameter by the 
relation Ha = y/ NRe which gives the typical value Ha ~ 700 for the present experiment. 
The influence of the electrical boundary conditions at the end faces is determined by this 
number and the conductivity ratio K = where a w is the electrical conductivity of 

the walls and l w their thickness. A very detailed numerical study on this issue has been 
performed in the laminar situation on a flow which geometry is very similar to the present 
one in [ijj]. As far as turbulent flows are concerned, Eckert et al. stress the importance of 
the product KHa, which measures the fraction of the electrical current that leaks from the 
Hartmann layers into the electrically conducting walls. With 12 mm thick end faces K is 
about 0.05 so that KHa ~ 35: the walls have to be considered as electrically conducting. 
These boundary effects are of crucial importance in experiments where the boundary layers 
control the rate of turbulence of the flow. However, in the present experimental setup the 
fluid is forced inertially and we have checked that the magnetic field has only little influence 
on the mean flow (see section Hi Bl below). Furthermore, the high value of the Hartmann 
number implies that the dominant balance in the bulk of the flow is not between the Lorentz 
and viscous forces but between the Lorentz force and the inertial term of the Navier-stokes 
equation: in the following, it is the interaction parameter and not the Hartmann number 
that leads to a good collapse of the data onto a single curve. For these reasons the present 
experimental device allows to study the influence of a strong magnetic field on the forced 
turbulence generated in the central shear layer, avoiding any boundary effect. 

C. Probes and measurements 

The propellers are driven at constant frequency. The current provided to the DC motors 
is directly proportional to the torques they are applying. It is recorded to access the 
mechanical power injected in the fluid. 

As liquid metals are opaque, the usual velocimetry techniques such as Laser Doppler 
Velocimetry or Particle Image velocimetry cannot be used. Hot wire anemometry is difficult 
to implement in liquid metals even if it has been used in the past [2|. Other velocimetry 
techniques have been developed specifically for liquid metals. Among those are the potential 
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FIG. 2: Schematics of the potential probes. These probes are vertical and the electrodes are 
positioned in the mid-plane of the tank, where the shear layer induces strong turbulence. 

probes. They rely on the measurement of electric potential differences induced by the motion 

nnn 

of the conducting fluid in a magnetic field pj, H3J, 114(1 • The latter can be applied locally with 
a small magnet or at larger scale as in our case. We built such probes with 4 electrodes 
(Fig. [2]). The electrodes are made of copper wire, 1 mm in diameter, and insulated by a 
varnish layer except at their very tip. The electrodes are distant of I ~ 3 mm. The signal 
is first amplified by a factor 1000 with a transformer model 1900 from Princeton Applied 
Research. It is further amplified by a Stanford Research low noise preamplifier model SR560 
and then recorded by a National Instrument DAQ. Note that because of the transformer the 
average potential cannot be accessed, so that our study focuses on the turbulent fluctuations 
of the velocity field. 

For a steady flow, and assuming j = 0, Ohm's law gives V0 = v x B so that the electric 
potential difference between the electrodes is directly related to the local velocity of the 
fluid. One gets 5<p = j v x B ■ dl integrating between the two electrodes. If B is uniform 
then 5<p = v±Bl where v± is the component of the velocity orthogonal to both the magnetic 
field and the electrodes separation. In the general time- dependent case, the link is not so 
direct. Using Coulomb's gauge and taking the divergence of Ohm's law, one gets : 



where io = V x v is the local vorticity of the flow [14|. The measured voltage depends on 
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FIG. 3: (a) Time series from a potential probe for f ro t = 20 Hz and an applied magnetic field 
of 178 G. (b) Corresponding power spectrum density. The dashed line is a /~ 5 / 3 scaling and the 
mixed line a / _11//3 scaling. These straight lines are drawn as eye guides. 

the vorticity component parallel to the applied magnetic field, i.e. to gradients of the two 
components of velocity perpendicular to B . The relation between the flow and the potential 
is not straightforward but the potential difference can be seen as a linearly filtered measure- 
ment of the velocity fluctuations. For length scales larger than the electrode separation, the 
potential difference can be approximated by the potential gradient, which has the dimension 
of vB . The spectral scaling of V<j)/B is expected to be the same as that of the velocity, 
i.e. the k~ 5 ^ 3 Kolmogorov scaling. For smaller scales, some filtering results from the finite 
size of the probe. For scales smaller than the separation I, the values of the potential on the 
two electrodes are likely to be uncorrelated: if these scales are also in the inertial range, one 
expects the spectrum of the potential difference to scale as the potential itself, i.e. fc -11 / 3 
due to the extra spatial derivative. Assuming sweeping of the turbulent fluctuations by the 
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average flow or the energy containing eddies [15] and a k~ 5 ^ 3 Kolmogorov scaling for the 
velocity spectrum, we expect the temporal spectrum of the measured potential difference to 
decay as f~ 5 ^ 3 for intermediate frequencies and as f~ n ^ 3 for high frequencies in the inertial 
range. 

An example of a measured time series of the potential difference is shown in figure [3] 
together with its power spectrum density. For this dataset, one expects a signal of the order 
of 2iTRf rot B l ~ 20 «V which is the right order of magnitude. Because of the quite large 
separation of the electrodes, the cutoff frequency between the f~ 5 ^ 3 and the f" 11 ^ 3 behaviors 
is low so that no real scaling is observed but rather trends. At the highest frequencies, the 
decay is faster than f~ n ^ 3 . 



Bolonov et al. [16( had a rather empirical approach to take into account the filtering 
from the probe. Assuming that the velocity spectrum should decay as f~ 5 ^ 3 they observed a 
spectral response of the potential probe which displays an exponential decay exp(— ///0.6m). 
The factor 0.6 is most likely dependent on the geometry. We reproduced their analysis in 
figure HI In the inset is displayed the spectrum of figure [3] multiplied by the expected 
j-5/3 sca li n g. The decay is seen to be exponential from about 25 Hz to 300 Hz (at higher 
frequencies the signal does not overcome the noise). The characteristic frequency of the 
decay can be extracted and plotted as a function of the rotation frequency of the propeller 
(the local velocity is expected to scale as Rf ro t)- A clear linear dependence is observed, in 
agreement with the results of Bolonov et al. The cutoff frequency is about twice f rot . It 
is not very high due to the rather large size of the probe. Nevertheless the most energetic 
length scales are resolved in our measurement. In the following we show only direct spectra, 
and no correction of the filtering is attempted. 

We can conclude that although some filtering is involved, the measurement of potential 
differences gives an image of the spectral properties of the velocity fluctuations. Any change 
in the spectral properties of the flow in the vicinity of the probe will thus be visible on the 
spectrum of the potential difference. 

The potential probes are quite large in order to fit a gaussmeter Hall probe in the vicinity 
of the electrodes (see Fig. [2]). These probes are connected to an F.W.Bell Gaussmeter model 
7030 that allows to measure the induced magnetic field down to a few tenths of Gauss. The 
proximity between the velocity and magnetic field measurements allows to study the possible 
correlations between these two fields. 
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FIG. 4: (a) Filtering effect due to the probe geometry. - Inset: semilog plot of the power spectrum 
density (PSD) of the potential difference compensated by / 5//3 . The dashed line is an exponential 
fit. - Main figure: lower curve, PSD; upper solid line: PSD corrected from the exponential decay 
fitted in the inset. Dashed line: /~ 5 ^ 3 decay. Same data set as that of the previous figure, (b) 
Cutoff frequency of the probe response for the probe used in (a) and for various rotation frequencies 
of the propeller. The dashed line is a linear fit. The applied magnetic field is 178 G. 



II. EFFECT OF THE APPLIED MAGNETIC FIELD ON THE TURBULENCE 
LEVEL 



A. Velocity field 

We first focus on the fluctuation level of the velocity field accessed through measurements 
of potential differences. The evolution of the rms value of the potential difference is displayed 
in figure [5] as a function of the applied magnetic field and of the rotation rate. From equation 
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FIG. 5: Evolution of the rms value of the potential difference as a function of the applied vertical 
magnetic field, for different values of the rotation frequency. •: f ro t = 5 Hz, A: 10 Hz, □: 15 Hz, 
o: 20 Hz and *: 30 Hz. The solid lines correspond to the azimuthal potential difference and the 
dashed line to the radial potential difference for the same four electrode potential probe. 

([2]), the potential difference should behave as S(fi rms oc B v/l. For low values of the applied 
magnetic field, the interaction parameter N is low: the magnetic field has almost no effect 
on the flow, and the velocity scales as Rf ro t- Thus, for high f rot and low B (i.e. low N), 
5 (fi rms sould be linear in both f rot and B . The upper curve corresponds to the highest 
velocity f ro t = 30 Hz. For low Bo there is a linear increase of 5<fi rms . Then it seems to 
saturate. For smaller rotation rates, the linear part gets smaller and the saturation region 
gets wider. Eventually, for f rot = 5 Hz, the potential decays for the highest values of the 
applied magnetic field. This demonstrates that there is a strong interaction between the 
magnetic field and the flow. 

To investigate in more details the scaling properties of the potential difference, 5(fi rms / f rot 
is displayed versus B Q in figure El One can clearly see that as f rot decreases, the potential 
deviates from the linear trend for lower and lower values of the magnetic field. 

Figure [7] shows 6(fi rms /B Q as a function of f rot . Here a linear trend is observed for large 
f rot . As the external magnetic field is increased, the fluctuations of the potential are damped 
and the linear trend is recovered for increasingly high values of the rotation rate. 

All that information can be synthesized by plotting the dimensionless potential 
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FIG. 6: Evolution of the rms value of the azimuthal potential difference normalized by the rotation 
frequency, as a function of the applied vertical magnetic field. •: f ro t = 5 Hz, A: 10 Hz, □: 15 Hz, 
o: 20 Hz and *: 30 Hz. The dashed line is a linear trend fitted on the first four points of the 30 
Hz data. 

, x 1(T 3 



0.8 
m° 0.6 

CO 

E 

<t~ 0.4 
0.2 



0. 








10 15 20 25 30 
f , [Hz] 
rot 1 ' 



FIG. 7: Evolution of the rms value of the azimuthal potential difference normalized by the imposed 
vertical magnetic field, as a function of the rotation frequency of the propeller. •: Bq = 356 G, 
A: 712 G, □: 1070 G, o: 1420 G. The parallel dashed lines are used as eye guides. The upper one 
corresponds to a linear fit of the data at Bq = 356 G. 



This quantity can be understood as the velocity fluctuations of the flow normalized by the 
forcing velocity of the propeller. This quantity is displayed as a function of the interaction 



parameter N 



oLBl 



in figure [8j In this representation, the data collapses fairly well on 



2npRf ro t 

a single master curve for high rotation rates. The damping of the turbulent fluctuations 
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FIG. 8: Evolution of the dimensionless potential a* (see text) as a function of the interaction 
parameter N. (a) and (b) correspond to the azimuthal potential difference taken from two different 
datasets. The potential probes are similar but not perfectly identical. The rotation frequencies 
reach the highest values in (a) and the lowest in (b). The representation is semilogarithmic. 

can reach one order of magnitude for N close to 1. In fig. [H(b) - that corresponds to lower 
rotation rates - there is a slight and systematic drift of the curves with the rotation rate. 
This indicates a slight dependence in Reynolds number that comes likely from the fact 
that the flow is not fully similar with f rot for low values of this rotation rate. The main 
dependence is clearly in the interaction parameter. The velocity fluctuations are seen to 
decay exponentially with N for values of N up to order 1. For higher N the data is affected 
by the noise. The velocity fluctuations are so strongly damped that the signal to noise ratio 
decreases significantly, as can be seen on the spectra in the following sections. The decay rate 
is about 2.5 in fig. [8](a) and 3.5 in (b). The difference may come from geometrical factors of 
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the probes which are not exactly similar in both datasets, or from a slight mismatch between 
the positions of the probes in the two datasets. We observed a similar collapse of potential 
data with N in a different experiment jl7|. This older experiment was smaller in size, with 
only one propeller and smaller magnetic field. Only the beginning of the exponential decay 
could be observed in that case. 



B. Induced magnetic field 




FIG. 9: (a) Evolution of the rms value of the induced radial magnetic field as a function of the 
applied vertical magnetic field, for different values of the rotation frequency. •: f ro t = 5 Hz, A: 
10 Hz, □: 15 Hz, o: 20 Hz and *: 30 Hz. (b) Evolution of 

brrms/ BqRtti as a function of the 
interaction parameter N. The representation is semilogarithmic. 



The induced magnetic field is of the order of one percent of the applied field. This low 
value is due to the low magnetic Reynolds number which is at best of order 1. For low Rm 
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the induction equation reduces to 



B ■ Vv + r]Ah = (9) 

(assuming that Bo is uniform). The induced magnetic field thus reflects the velocity gradients 
in the direction of the applied magnetic field. We use it as a second tool to investigate the 
statistical properties of the flow. 

The fluctuation level of the radial magnetic field b rrms is displayed in figure [9](a). Its 
evolution with Bq and f ro t is strongly similar to that of the potential differences, as expected 
from the previous arguments. The azimuthal component displays the same behavior (not 
shown) . 

For low values of the applied magnetic field, the flow is not affected very much by the 
Lorentz force and from equation one expects the induced magnetic field amplitude to 
scale as: b oc B^Rm (see [3] for example). From the previous section, the amplitude of 
the velocity field decays exponentially with N. We thus expect b rms /BoRm to have the 
same qualitative behavior. This is what is observed in figure Mjo). The various datasets 
are seen to collapse on a single master curve for N up to 0.5. For higher N, the sensitivity 
of our gaussmeter is not high enough for the signal to overcome the electronic noise. As 
for the velocity, the curve for the lowest f rot is below all the others, which confirms the 
slight dependence with Rm observed on the potential measurements. At a given value of 
N, the collapse of the measurements shows that the fluctuations of the induced magnetic 
field are indeed linear both in applied magnetic field and in magnetic Reynolds number: 
b rms ~ 0.1B Rm. 

We have also measured the average value of the induced magnetic field. From equation 
([§]) it is linked to the vertical gradients of the time-averaged velocity field. (b r ) / B^Rm 
and {bg) / BqRtji are shown in figure [10] as a function of N. For both components, this 
representation leads to a good collapse of the datasets. At a given value of N, the collapse 
shows again that the average induced magnetic field is linear in Rm and B with ~ 
0.01B Rm for the data shown here. In an ideal Von Karman experiment, when the two 
propellers counter-rotate at the same speed, the time averaged flow is invariant to a rotation 
of angle tt around a radial unit vector taken in the equatorial plane (denoted as e r on figure 
[Q. If the applied field were perfectly symmetric and the probe positioned exactly in the 
equatorial plane, this symmetry should lead to (b r ) = 0. However, the introduction of the 

16 



0.015 
0.01 

E 

c 0.005 

O 

CO 



A 
V 





-0.005 



-0.01 



• 


5 Hz 


▲ 


10 Hz 


■ 


15 Hz 




20 Hz 


★ 


30 Hz 



(a) 



0.2 0.4 0.6 O.i 
N 



0.01 
0.01 
QC E 0.008 

o 
CD 

~ 0.006 



"v 0.004 



0.002 




• 


5 Hz 


▲ 


10 Hz 


■ 


15 Hz 


♦ 


20 Hz 


★ 


30 Hz 



(b) 



FIG. 10: (a) Evolution of 
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as a function of the interaction parameter N. (b) same for 
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probe breaks the symmetry and neither the mecanical device nor the applied magnetic field 
are perfectly symmetric. It has been observed in our setup and in the VKS experiment 
that the measurement of (b r ) is extremely sensitive to the position of the probe (private 
communication from the VKS collaboration). For example, a slight mismatch of the propeller 
frequencies or of the relative positions of the probe and mid-plane shear layer can lead to 
strong changes in the mean value of the magnetic field. This may be the reason why (b r ) 
is not zero here. Nevertheless there is a systematic change of the average with N which is 
consistent across the different values of the velocity and of the applied magnetic field. It 
indicates a change of the time- averaged flow in the vicinity of the probe. To the shear layer 
lying in the equatorial plane corresponds strong radial vorticity, orthogonal to the applied 
magnetic field. The applied strong magnetic field will impact the shear layer, as it tends to 
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elongate the flow structures along its axis. Even a small change in the shear layer geometry 
affects strongly the measured (b r ). Here we see that its sign is changed at high N. 



The time-averaged induced azimutha! 



differential rotation of the propellers 
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magnetic field (bo) is due to u-eSect from the 



181 ] . It is seen to decay slightly (about 35%) with 



N. This effect may be due to some magnetic braking that leads to an elongation of the shear 
layer and thus to weaker differential rotation in the vicinity of the mid-plane. Nevertheless 
this effect is limited and we expect the average large scale structure of the flow to remain 
almost unchanged. The small change of the large scale flow is not the reason for the strong 
damping of the turbulent fluctuations by an order of magnitude. 



C. Injected mechanical power 



300 




FIG. 11: Injected mechanical power as a function of the rotation frequency of 
the propellers. The circles correspond to measurements without magnetic field, and 
the triangles correspond to measurements at constant rotation frequency and -Bo = 
90, 180, 530, 705, 880, 1000, 1230, 1240, 1320, 1410, and 1500 Gauss. The solid line has equation 
y = 0.032/^ which comes from the turbulent scaling law. The inset is a log-log representation. 

An interesting issue in MHD turbulence is to understand how the injected mechanical 
power is shared between ohmic and viscous dissipations. The torques T\ and T2 provided 
by the motors can be accessed through measurements of the current delivered to them. The 
injected mechanical power is then e = (T x + T 2 )2iif rot . This quantity has been measured as a 
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function of the rotation frequency with and without applied magnetic field. The results are 
drawn on figure [TTJ without magnetic field, the injected power follows the turbulent scaling 
law e ~ fr 0t . More surprisingly, we notice that at a given rotation frequency the experimental 
points corresponding to different amplitudes of the magnetic field are indistinguishable. This 
observation confirms that no change of the global structure of the flow is induced by the 
magnetic field. The fact that the injected power is independent of the applied magnetic 
field seems to be in contradiction with results from numerical simulations where an ohmic 
dissipation of the same order of magnitude as the viscous one is observed when a strong 
magnetic field is applied (the ohmic dissipation is around three times the viscous dissipation 
for iV = 1 in Burattini et al. [8|). One could argue that there may be a large ohmic dissipation 
compensated by a drop in the viscous dissipation: the energy flux in the turbulent cascade 
would then be dissipated mainly through ohmic effect without changing the overall injected 
power. However, a rough estimate of the ohmic dissipation Dj gives values which are rather 
low: with the maximum value of the induced magnetic field b ~ O.lBoRm, and assuming 
that the magnetic field is dissipated mostly at large scale, one gets: 

c ,i^«, w (10 , 

where we used the values Rm = 1, and B = 1500G. This estimate is much lower than 
the injected power and we may expect ohmic dissipation to remain negligible compared 
to viscous dissipation, even at order one interaction parameter. However, one also needs 
to know the current that leaks through the boundaries to evaluate the additional ohmic 
dissipation that takes place inside the walls. As these effects are difficult to quantify, we are 
not able to determine the ratio of ohmic to viscous dissipation. Nevertheless it is interesting 
to notice that the injected mechanical power remains the same when a strong magnetic field 
is applied, although velocity fluctuations are decreased by a factor 10 in the mid-plane of 
the tank. 



D. Development of anisotropy 



The application of a uniform magnetic field on a turbulent flow is known to elongate 
the flow structures in the direction of the applied field [3(. In decaying turbulent flows, 
this effect eventually leads to the bidimensionalization of the flow [19|. As far as forced 
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turbulence is concerned, only numerical simulations have demonstrated this effect [7]. The 
measurement of both the induced magnetic field and the electric potential allows to quantify 
the elongation of the turbulent structures in the z direction: on the one hand, equation (J7J) 
links the electric potential to the vertical component of vorticity, i.e. to horizontal gradients 
of velocity. On the other hand the induced magnetic field is related to vertical gradients of 
velocity through equation ([9]). We expect then that 



Ab 1 d\\v 



(11) 



Acj) j) d±v 

where du and d± denote derivatives in directions parallel and perpendicular to the applied 
magnetic field. As we cannot access experimentally the order of magnitude of the Laplacians, 
we define the quantity a = r\ b s r £ ms which we expect to give a crude estimate of the ratio of 
the vertical to the horizontal gradients of velocity. It is somewhat related to the parameter 
G\ defined in Vorobev et al. [7J. This anisotropy parameter a is represented as a function 
of the interaction parameter N in fig. [12] for different values of the rotation frequency: it 
decreases from about 4 until it saturates around 1.5. This decrease of the parameter a by a 
factor about 3 when a strong magnetic field is applied provides evidence for the elongation 
of the flow structures in the z direction: the derivatives of the velocity field in the direction 
of B are much smaller than its derivatives in directions perpendicular to B . Although the 
turbulence becomes more anisotropic, it remains three dimensional even for the highest value 
of the interaction parameter reached in this experiment. This is due to the fact that the 
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forcing imposed by the propellers is three dimensional and prevents the flow from becoming 
purely two dimensional. 

III. TEMPORAL DYNAMICS 

We observe that the turbulent fluctuations are being damped by magnetic braking when 
a strong magnetic field is applied to homogeneous and nearly isotropic turbulence. An 
interesting question is to know how this damping is shared among scales. In this section we 
study the evolution of the power spectrum densities of the potential difference and induced 
magnetic field. 

A. Potential 

We show in figure [TBI the power spectrum densities (PSD) of the dimensionless potential 
v* = <p/ B lRf rot . For a rotation frequency of 20 Hz, the spectra are seen to decay as the 
magnetic field is increased, the shape of the different spectra remaining the same (fig. ITBT c)). 
The decay of the PSD is at most of a factor 10. This dataset reaches a maximum value of the 
interaction parameter about 0.4. For the lowest displayed value of the rotation frequency 
frot = 5 Hz (fig. [131(a)). N reaches 1.5. As the interaction parameter increases, two distinct 
regimes are observed: first an overall decay of the PSD and second a change in the shape 
of the PSD. The highest frequencies are overdamped compared to the lowest ones: for the 
highest value of N, the PSD decays by about 6 orders of magnitude at twice the rotation 
frequency, whereas it decays by only two orders of magnitude at low frequency. We already 
observed the first regime in a previous experiment performed on a different flow jl7| . The 
interaction parameter was not high enough in this experiment to observe the second regime. 

To quantify more precisely the relative decay, we plot in figure [14] the ratio of the PSD 
of v* over the PSD at B = 178 G (and at the same rotation rate of the propeller). At 
the smallest values of N, the ratio weakly changes across the frequencies but it decays with 
N. When N is increased over approximately 0.1 an overdamping is observed at the highest 
frequencies. This extra damping can be qualitatively characterized by a cutoff frequency, 
which decreases extremely rapidly and seems to reach the rotation frequency for N ~ 0.3. 
Above the cutoff frequency, the decay rate seems to behave as a power law of the frequency. 
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for the radial potential difference. The subfigures correspond to different rotation rates of the 
propellers: (a) f ro t = 5 Hz, (b) 10 Hz, (c) 15 Hz. In each subfigure, the different curves correspond 
to the different values of the applied magnetic field B = 178, 356, 534, 712, 890, 1070, 1250, 1420 
and 1600 G. They are naturally ordered from top to bottom as the magnetic field is increased. The 
noise part of the spectra has been removed to improve the clarity of the figures. 

The exponent of this power law gets more and more negative as N increases and seems to 
reach —5 at the highest value of N displayed here (f rot = 5 Hz, B = 1600 G, N = 1.5). 

In figure [151 we gathered the dimensionless spectra at various rotation rates that corre- 
spond to the same interval of N. The spectra are collapsing fairly well onto each other. A 
little bit of scatter is observed, most likely due to the slight dependence on Rm described 
previously. The shape of the spectra is essentially a function of the interaction parameter. 

In experiments on the influence of a magnetic field on decaying turbulence, the velocity 
spectrum goes from an f~ 5 ^ 3 to an f~ 3 behavior as the interaction parameter is increased. 
This —3 exponent is attributed either to two-dimensional turbulence or to a quasi-steady 
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FIG. 14: Evolution of the shape of the dimensionless potential's power spectrum. The spectra of 
figure [13] are divided by the spectrum obtained at the smallest value of the applied magnetic field 
Bq = 178 G and at the same f ro t- The subfigures correspond to different rotation rates of the 
propellers: (a) f ro t = 5 Hz, (b) 10 Hz, (c) 15 Hz. 

equilibrium between velocity transfer and ohmic dissipation. As far as forced turbulence 
is concerned, we observe a strong steepening of the velocity spectrum, with slopes already 
much steeper than f~ 3 for N — 1. We do not observe any signature of this quasi-steady 
equilibrium or of 2D turbulence. Once again, this comes from the three-dimensional forcing 
of the propellers which prevents the flow from becoming purely 2D. 

Using the full ensemble of datasets at the various rotation rates one can interpolate the 
evolution of the shape of the dimensionless spectra as a function of both / / f rot and N. The 
result is shown in figure [161 This representation summarizes all previous observations. As 
iV increases, first there is a self similar decay of the spectra. When N reaches approximately 
0.1, an additional specific damping of the high frequencies is observed. The high frequency 
spectrum gets extremely steep. This extremely steep regime covers the full inertial range 
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FIG. 15: Evolution of the shape of the normalized potential's power spectrum as a function of the 
interaction parameter. Colors corresponds to the different rotation speeds: black 3 Hz, blue 5 Hz, 
green 7.5 Hz and red 10 Hz. Each subfigure corresponds to data restricted to the specified interval 
of N. The dashed line corresponds to the spectrum at the smallest non zero value of N. 




FIG. 16: Evolution of the shape of the power spectrum of the normalized potential as a function 
of the interaction parameter. Data for 5, 10, 15 and 20 Hz have been used for this representation. 
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for N close to 1. 

B. Induced magnetic field 




FIG. 17: Evolution of the power spectrum of the dimensionless induced magnetic field . 

BoRm 

The subfigures correspond to different rotation rates of the propellers: (a) 10 Hz, (b) 15 Hz and (c) 
20 Hz. In each subfigure, the various curves correspond to different applied magnetic fields. The 
curves are naturally ordered from top to bottom as the magnetic field is increased. The results 
for the following values of B are displayed: B = 178, 356, 534, 712, 890, 1070, 1250, 1420 and 
1600 G. The noise part of the spectra has been removed for the clarity of the figures. 

The same analysis can be performed on the induced magnetic field. The dimensionless 
spectra of b r are shown in figure [T71 One concern is that the induced magnetic field is low, so 
that the signal to noise ratio of the gaussmeter is not as good as that of the potential probes. 
The magnetic field spectrum reaches the noise level at a frequency which is approximately 
2f rot . This is also related to the fact that the scaling law expected from a Kolmogorov-like 
analysis for the magnetic field is much steeper than the one of the velocity (/~ n / 3 in the 
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dissipative range of the magnetic field, at frequencies corresponding to the inertial range of 
a velocity scaling as / _5//3 ). Nevertheless, the same two regimes are observed: at low N 
the spectra are damped in a self similar way. At large N the high frequencies seem to be 
overdamped. However, the picture is not so clear because of the low signal to noise ratio. 
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FIG. 18: Evolution of the shape of the power spectrum of the induced magnetic field b r . The 
spectra of figure [T7] are divided by the spectrum obtained at the smallest value of the applied 
magnetic field Bq = 178 G and at the same f ro t- The subfigures correspond to different rotation 
rates of the propellers: (a) 10 Hz, (b) 15 Hz and (c) 20 Hz. 



This last point is better seen on the ratio of the spectra in figure [181 The ratio corre- 
sponding to the highest applied magnetic field is decaying at large frequencies. 



C. Coherence between velocity and magnetic field 

The experiment has been designed to record the potential difference and the magnetic 
field in the vicinity of the same point. Figure [19] displays the spectral coherence between 
these two quantities. In this figure we have identified the azimuthal potential difference with 
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(a) f[Hz] (b) '[Hz] 

FIG. 19: Spectral coherence between the velocity and the induced magnetic field, (a) coherence 
between b r and v r or vg. The inset recalls the spectra of b r and vg. (b) coherence between bg and 
v r or vg. The inset recalls the spectra of bg and v r . The data correspond to f ro t = 15 Hz and 
B = 356 G. 

the radial component of velocity v r and the radial potential difference with the azimuthal 
component of velocity Vg, even though this identification may be somewhat abusive. We 
recall that the coherence is one when both signals are fully correlated at a given frequency 
and zero if they are uncorrelated at this frequency. We have plotted in the insets the 
power spectrum of each signal. The magnetic field spectrum falls below the noise level at 
about 2 f ro t, which explains why all coherence curves go to zero above 2f rot (there is no 
magnetic signal at these frequencies). For lower frequencies, a small but non-zero value of 
the coherence is observed for the following pairs: (b r ,v r ) with a coherence level above 0.1, 
(b r , vg) with a coherence close to 0.05 and (bg, vg) with a coherence level barely reaching 0.1. 
The coherence value changes a little for other values of f rot , but the global picture remains 
the same, with (b r ,v r ) being the most coherent. 

The time correlations are shown in figure [201 No clear correlation is observed except for 
the pair (b r ,v r ) which displays a small peak reaching 0.1. The weak coherence observed in 
the previous figure is barely seen here, probably because of an insufficient convergence of 
the correlation function. 

At the low values of magnetic Reynolds number attained in our experiment, the magnetic 
field is diffused through Joule effect, so that its structure is mainly large scale. The Reynolds 
number being large, the velocity fluctuations develop down to much smaller scales. Because 
of the strong ohmic diffusion, the magnetic field is expected to be sensitive mostly to large- 



27 



0.1 

O 

" a, 

:> 

O 

-0.1 

-5 5 

time (s) 

v e 

~-5~~ 0~ 5 

time (s) 

FIG. 20: Correlation coefficients between the velocity and the induced magnetic field, (a) Correla- 
tion between b r and v r (upper curve) or vg (lower curve), (b) Correlation between bg and v r (lower 
curve) or v r (upper curve). The data correspond to f ro t = 15 Hz and Bq = 356 G. 

scale and low frequency fluctuations of the velocity field. For example one could expect 
to observe a correlation or coherence between vg and bg: fluctuations of the differential 
rotation can induce fluctuations in the conversion of the axial Bo into azimuthal magnetic 
field through uo effect. As far as b r is concerned, the poloidal recirculation bends the vertical 
field lines towards the exterior of the tank in the vicinity of the mid-plane, a process which 
induces radial magnetic field from the vertical applied field Bo- Fluctuations of this poloidal 
recirculation thus directly impacts the radial induced magnetic field, hence the correlation 
between v r and b r . Although such low frequency coherence is indeed observed in figure [19] 
its amplitude remains very low, so that there is almost no correlation between the velocity 
field and the magnetic field measured in the vicinity of the same point. 

IV. DISCUSSION AND CONCLUSION 

The effect of a strong magnetic field on forced turbulence is studied experimentally with 
potential probes and induced magnetic field mesurements. The velocity fluctuations are 
strongly damped as the applied magnetic field increases: for N ~ 0.6, the turbulence inten- 
sity in the mid-plane of the tank is decreased by an order of magnitude. As a consequence, 
the standard deviation of the induced magnetic field - normalized by BqRiti - also diminishes 
by a factor ten. The spectrum of the non-dimensional potential v* is affected in two differ- 
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ent ways by the magnetic field: for low values of the interaction parameter the spectrum 
decays uniformly at all frequencies, its shape remaining the same. For higher values of N, 
we observe an overdamping of the high frequencies. The same effect is seen on the induced 
magnetic field spectra. 

We identify several features which highlight the very different behaviors of forced and 
decaying turbulence when they are subject to a strong magnetic field: decaying turbulence 
is thought to evolve towards a bidimensional structure. Its velocity spectrum displays a —3 
exponent which can be attributed either to this bidimensionalization or to a quasi-steady 
equilibrium between velocity transfer and ohmic dissipation. In the present experiment, 
no such —3 exponent is observed, and the spectra are much steeper (exponent —5 to —6) 
for values of N higher than 0.5. Moreover, we have introduced a parameter a = ri ^ ms to 

0<Prms 

quantify the anisotropy of the turbulence in the mid-plane of the tank. When N increases, 
the decrease of this parameter by a factor 3 is the signature of the elongation of the flow 
structures along the applied magnetic field. However, the flow always remains 3D since a is 
non-zero even for the highest value of the interaction parameter reached in this experiment. 
These differences between forced and decaying turbulence come from the 3D forcing imposed 
by the propellers, which rules out the possibility of a 2D statistically steady state of the flow. 
It would be interesting to perform the same kind of experimental study of the anisotropy 
with other forcing mechanisms, such as current-driven MHD flows or turbulent thermal 
convection in a liquid metal. 

We have studied the evolution of the injected mechanical power as the applied magnetic 
field increases and found almost no influence of the latter: the injected mechanical power 
remains the same although velocity fluctuations are decreased by an order of magnitude in 
the central shear layer. 

Finally, we stress the poor level of correlation between the velocity and induced magnetic 
fields measured in the vicinity of the same point, and attribute it to the scale separation 
between the two fields. 

The strong damping of turbulent fluctuations by the magnetic field can be invoked as 
a saturation mechanism for turbulent dynamos: in a dynamo experiment, one observes 
spontaneous generation of magnetic field when the magnetic Reynolds number is above a 
critical value Rm c . If the turbulent fluctuations are involved in the generating process of the 
magnetic field (through au or a 2 mechanisms for instance), there is a critical level of rms 
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turbulent fluctuations a vc above which magnetic field is generated (a v = a vc for Rm = Rm c ) . 
For Rm > Rm c , the initial level of turbulent fluctuations is above a vc , and the magnetic field 
grows exponentially from a small perturbation: the interaction parameter increases, and the 
turbulent fluctuations are damped according to figure El An equilibrium is reached when the 
damping is such that the rms turbulent fluctuations are reduced to a vc . For small values of 
N we observed that a v (N) = a v {N = 0)e^ N ~ a v (N = 0){1-^N), with 2.5 < 7 < 3.5. The 
saturated value N sat of N follows from the equality cr v (N sat ) = <r vc — a v (N — 0)(1 — -yN sat ). 
As cr v is proportional to Rm, we get N sat = Rm ^^ nc , hence the following scaling law for the 
magnetic field: 



2 2npRf rot Rm — Rm c 



1 - : L 7 Rm ; (i2) 



This is the turbulent scaling law for the saturation of a dynamo, which was originally 
described by Petrelis et al. (see [2^] for instance). Although the exact geometry of the 
experimental setup and large scale magnetic fields are different from that of the present 
experimental study, results from the VKS dynamo can be used to test this relationship: 
using R = 15.5 cm, L = 21 cm, f rot = 16 Hz, p = 930 kgm~ 3 , a = 9.5 10 6 Vt^m' 1 and 
7 = 3, the computed magnetic field amplitude is B sat ~ 290 G for Rm f [ r ^ mc = |- This is the 
right order of magnitude: the amplitude of the VKS dynamo field measured in the vicinity 



2l|: figure 3(b)). 



of the axis of the cylinder is approximately 150G for this value of Rm 
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